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Abstract 

We consider evolutionary models for the population of short-period (Por 
(LMBHB) and compare them with observations of soft X-ray transients (SXT). Evolution of LMBHB is determined 
by nuclear evolution of the donors and/or orbital angular momentum loss due to magnetic braking by the stellar 
wind of the donors and gravitational wave radiation. We show that the absence of observed stable luminous LMBHB 
implies that upon RLOF by the low-mass donor angular momentum losses are substantially reduced with respect 
to the Verbunt & Zwaan "standard" prescription for magnetic braking. Under this assumption masses and effective 
temperatures of the model secondaries of LMBHB are in a satisfactory agreement with the masses and effective 
temperatures (as inferred from their spectra) of the observed donors in LMBHB. Theoretical mass-transfer rates 
in SXTs are consistent with the observed ones if one assumes that accretion discs in these systems are truncated 
("leaky"). We find that the population of short-period SXT is formed mainly by systems which had unevolved or 
slightly evolved {Xc ^ 0.35) donors at the Roche-lobe overflow. Longer period (Porb — (0.5 — 1) day) SXT might 
descend from systems with initial donor mass about 1 M© and Xc ^ 0.35. It is unnecessary to invoke donors with 
almost hydrogen-depleted cores to explain the origin of LMBHB. Our models suggest that a very high efficiency of 
common-envelopes ejection is necessary to form LMBHB, unless currently commonly accepted empirical estimates of 
mass-loss rates by winds for pre-WR and WR-stars are significantly over-evaluated. 
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1. Introduction 

Out of twenty dynamically-confirmed black-hole 
candidate X-ray binaries ten have K/M spectral 
type secondaries and orbi t al per iods ^ 1 day 



20061 ). All these sys- 



(jRemillard and McClintockl . 
tems are transient in X-ray. Their quiescent X-ray 

luminosity is ~ (310^° 310^^) ergs~^, while 

in the outburst luminosity may be by a factor 
~ 10^ — 10^ higher. These low-mass black-hole bi- 
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nary (LMBHB) X-ray systems belong to the class 
of "Soft X-ray transients" (SXT). Recent reviews 
of observational data o n SXT may be found e.g. in 
Charles and Coe (|2006h :l iRemillard and McClintockl 
(j2006r ). From the evolutionary stand-point SXT 
are considered as semidetached binaries in which 
matter is transferred to the compact object via ac- 
cretion disc. The commonly accepted model for the 
variability of SXT is based on the thermal -viscous 
insta bility of irradiated accretion discs (see [Lasota, 
I2OOIL and references therein). The estimated num- 
ber of S XT in the Galaxy ranges from several 
hundred (IChen et al.l . IiQQtI) to several thousand 
(Roman], 1998f ). Below, we will describe the models 
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Fig. 1. Initial — final mass relations for massive stars in bina- 
ries. The mass of stars at terminal main-sequence (TAMS), 
the initial mass of helium core Mcore.ii the pre-supernova 
mass Mpre_gN, and the black-hole mass Mbh a-re shown. 

of a population of short-period LMBHB, discuss 
their observational properties, and consider possi- 
ble existence of a population of faint LMBHB with 
quasi -stable, cold accretion discs ( Menou et al.l . 
19991). This review is ba sed mainly on the studies by 
Yungelson et al. ( 20061 ) (henceforth. Paper I) and 
Yungelson & Lasota (in prep.. Paper II). 



2. The model 

2.1. Formation of black holes with low-mass 
companions 

Modeling of the Galactic population of LMBHB 
includes two major steps: (a) modeling of the pop- 
ulation of zero-age binaries containing a black hole 
accompanied by a low-mass companion and (b) fol- 
lowing the evolution of the system to Hubble time. 

The population of zero-age LMBHB was com- 
puted (Paper I) using population synthesis code 
SEBA 



eccentricities of orbits S(e) = 2e. 

A summary of the assumptions on transforma- 
tions of stellar masses in the course of the evolu- 
tion is shown in Fig. [T] Black hole progenitors have 
Mq = 25 — 100 Mq. The algorithm for formation of 
black holes in the code follows the fall-back scenario 



(jFrver and Kalogeral 120011 ) under assumption that 
explosion energy is fixed at 10^*^ ergs, which is within 
the expected range but favours formation of rather 
massive black holes (up to ~ 15 M0). Nascent black 
holes r eceive an asymm etric kick at formation fol- 
lowing Paczvhskil (|l990ll velocity distribution with 
ay = 300 km s""'^, but scaled down with the ratio of 
the black hole mass to a neutron star mass. The pop- 
ulation of LMBHB is not sensitive to the assumed 
kick distribution, since scaled down with Afbh/-^^ns- 
ratio kick velocities are too small to disrupt close 
binaries in SN explosions. Test run for Maxwellian 
kick distribu t ion fo r pulsars with a = 265km/s 
(|Hobbs et al.l . l2005h confirmed this assumption. For 
the common- envelope phase we used the standar d 
prescription (jWebbinkl . Il984l : Ide Kool etali Il987t ). 
with the efficiency and structure parameters acE 
and A product acE x A = 2 (see discussion of the is- 
sue of acE X A in 35]) . For more detailed description 
of input parameters we refer the reader to Paper I. 

2.2. Angular momentum loss 

Since low-mass components in LMBHB have 
KV/MV spectral types, it is usually assumed that 
evolution of LMBHB is governed by angular mo- 
mentum losses (AML), like for cataclysmic vari- 
ables. In our computations we assumed that AML 
via magnetic stellar wind (MSW) follow s the model 
suggested bv lVerbunt and Zwaan (1981): 



din J 
dt 



X 10 



-28; 2 



(IPortegies Zwart and Verbunt . 1996 : Portegies Zwart and Yungelson 



19981 : iNelema ns et all . I2OOII : iNelemans et al.L l2004l ) 
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1 M^^^R. 



Q2/3 Mip ' 
(1) 

total mass 



with 250000 initial binaries with Mio > 25 Mq. 
Assumed time- and position-dependent Galactic 
star formation histor y is ba sed on the model of 
Boissier and Prantzod ( 1999f ). but for the inner 3 
kpc of the Galaxy we doubled the star formation 
rate given in the latter study to mimic Galactic 
bulge. We assumed a 50% binarity rate (2/3 of s tars 
in binaries), an IMF after iKroupa et al ( 19931 ). an 
initial distribution of separations in binaries (a) 
flat in logo between contact and 10® R©, a flat 
mass ratio distribution, and initial distribution of 



where Mi is mass of the primary, Mj 
of the system, P - orbital period, R - radius of the 
secondary, fc^ '--^ 0.1 - its gyration radius, / ~ 1 - 
a parameter derived from observations. Note, that 
this model extrapolates via almost two orders in 
magnitude s tellar rotation braking law derived by 
SkumanichI (|l972h for single field stars to rapidly 



rotating components of close binaries and assumes 
efficient spin-orbital coupling. 

Angular momentum loss via gravitational waves 
radiation (GWR) i s desc ribed by the standard 
iLandau and Lifshitz (1971) formula 
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Fig. 2. Fractions of low-mass bh+ms binaries formed at dif- 
ferent epochs in the hfetime of the Galaxy (soUd line) and the 
shape of star formation rate (SFR) adopted in present study 
(dashed line). At maximum, the SFR is ~ 15.6 Mq yr~^. 
Current SFR is ~ 1.9 M© yr"!. 



din J 
dt 



32 
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Ml M2 



(2) 



Below, we will refer to the model in which semide- 
tached systems obey Eqs. ([1]) and ([2]) as to a "stan- 
dard" model. 

2.3. Population of unevolved LMBHB 

With our assumptions we find that, within a Hub- 
ble time (13.5 Gyr), ~ 49 000 binaries that have or- 
bital periods below 2.0 day and contain black holes 
accompanied by main-sequence stars with M20 < 
1.6 M0 (for which magnetic braking is supposed to 
operate) were formed in the GalaxyUJ Out of them, 
~ 17000 were brought into contact by AML via 
MSW and/or GWR and ~ 14 000 evolved to shorter 
periods under "standard" assumptions on MSW (see 



Because black holes form in the first ~3 Myr af- 
ter their progenitor formation, formation history of 
bh+ms binaries strictly follows the star formation 
history (Fig. [J). Figure [3] shows relations between 
initial masses of components in model systems and 
their post-circularization orbital periocO]- Solid 
lines in the upper panel of the Figure outline the 



^ This number is ^^0.01% of all binaries with black-hole 
components formed in Hubble time. 

^ Every dot in the plot actually represents four or five sys- 
tems with similar combination of Mio, M20, Po, but born 
at different time. Since we present one random realization of 
the model, all numbers given are subject to Poisson noise. 
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Fig. 3. Upper panel: relation between orbital periods and 
masses of main-sequence companions of black holes after cir- 
cularization of the orbits. Systems with main-sequence com- 
ponents that may fill their Roche lobes in their main-se- 
quence lifetime or in the Galactic disc lifetime (13.5 Gyr) 
are located to the left of the solid curves that correspond to 
Mbh = 4M0 (left one) and 12 M© (right one), respectively. 
Lower panel: relation between initial masses of components 
in bh-f ms systems. 

mass and period range from which systems may be 
brought in the contact via action of AML, however, 
in "standard" model only systems that overflow 
Roche lobes at periods less than ~ 15 hr evolve to 
shorter periods (the rest have small helium cores 
and evolve to longer periods, forming at the end he- 
lium white dwarfs, if time span between contact and 
Hubble time is large enough; otherwise, at present 
their donors are subgiants). 

3. Evolution of LMBHB 

For simulation of the current Galactic popula- 
tion of semidetached LMBHB we convolved the 
above-described "underlying" population of bina- 
ries born at different epochs in the history of the 
Galaxy with the grids of evolutionary tracks for 
low-mass components in the binaries with differ- 
ent combinations of masses of components and 
post-circularization (initial) orbital periods. For 
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evolutionary computations, appro priately modified 
TWIN version (September 2003) of lEggletonI (|l97ll ) 
evolutionary code was used. Every "initial" system 
was "evolved" through 13.5 Gyr — Tq, where Tq is 
the moment of formation. In the Figures presented 
below, we show only systems that have at present 
q = M2/M1 > 0.02. At q < 0.02 the circularization 
radius of accretion stream becomes greater than the 
outer radius of the accretion disc, resonance phe- 
nomena in the disc become important and, though 
computations may be continued in the same fashion 
as before, until models cease to converge due to very 
low mass, actually it is unclear how mass transfer 
occurs (see Paper I for more details). The systems 
with q < 0.02 have Porb 2 hr and mass-transfer 
rates (< 10"^" Mq yr^^). In the "standard" model 
currently more than 75% of all low-mass ms+bh 
systems have q < 0.02 and have yet to be observed 
if mass-transfer occurs in them. 

3.1. Luminous persistent LMBHB? 



Figure H] shows the distribution of systems with 
q > 0.02 in the "standard" model over mass-transfer 
rates. The upper two panels show the break-down 
of the binaries that have stable or unstable accre- 
tion discs (370 and 2900 objects, respectively). Irra- 
diated discs are hot and stable if accretion rate (in 
M© yr-i) exceeds (|Dubus et al.l . ll999l) 



M, 



crit 



2.4 X 10-"AfB°*i?2jo 



C 



5 X 10- 



-0.5 



(3) 

where A/bh is the mass of the accretor in Mq, 
Rd,io is the disc outer radius in 10^° cm, and 
C is a measure of heating o f the disc by X-rays 
(jShakura and Sunvaevl 119731 ): for a given disc ra- 
dius R the irradiation temperature Ti„. is 



Mr/ 



(4) 



where cr is the Stefan- Boltzmann constan t. We use 
C = 5x10^"*, foUowing lDubus et al.l ( 19991 ) who have 
found that this value is consistent with properties 
of persistent low-mass X-ray sources, but actually C 



might e.g. vary in time l Esin et al. . 2000l ) 



However, no persistent LMBHB with Porb ^ 4 hr 
are observecPI. Stable model systems have an aver- 
age luminosity L ~ O.lA/c^ « 5 • 10^^ erg s^^ The 
accretion luminosity may be assumed to be emitted 
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In fact no LMBHB at all are observed below P„, 



4 hr. 



Fig. 4. Mass transfer rates vs. orbital period for the model 
population of semidetached low-mass binaries with assumed 
black hole accretors. Upper panel - stable systems in the 
model for a "standard" Verbunt-Zwaan assumption on AML 
via MSW. Middle panel - population of unstable systems in 
the same model. Lower panel - population for the model in 
which MSW does not operate in semidetached systems. In 
this model all binaries have unstable discs. The arrows mark 
approximate upper limits to the average mass-transfer rates 
for SXT under assumption of "leaky" discs, while plusses 
show estimates of rates based on assumed recurrence period 
of 30 years, except for V616 Mon and IL Lup whose recur- 
rence times are known to be ~ 50 and 10 years, respectively 
(see § 13.3.21 for details) . 



as that of a 1 keV blackbody, corresponding to the 
"high/soft" X-ray spectral state of LMBHB. Given 
their location in the Galaxy provided by the popula- 
tion synthesis code, we may roughly estimate inter- 
vening X-ray absorption. We found that the number 
of persistent systems expected above the complete- 
ness limit for the X-ray sources seen by the RXTE 
All Sky Monitor (about 3 mCrab or 0.2 counts/s) 
in the 2-10 keV band is 85. Setting the threshold 
ten times higher still yields 13 bright, persistent 
LMBHB with ASM fluxes comparable to e.g. SMC 
X-1. Noticeable number of stable "observed" persis- 
tent model sources remain in the model also under 
another assumptions on the spectrum of their emis- 
sion (see for details Paper I) . If the discs were not ir- 
radiated, no systems would be stable. However, both 
persistent neutron-star low-mass X-ray binaries and 
outbursting black hole transients show clear signa- 
tures of irradiation. Note that with parameter C = 



4 



5 X 10^^, only a 0.5% of the accretion luminosity 
(efficiency of 10%) is reprocessed in the disc. 

Reduction of dlnJ/dt by factor 2 compared to 
Eq. ^ does not change the results significantly (see 
Paper I) and points to the necessity of stronger re- 
duction of AML for semidetached systems. 



3.2. Reduced angular momentum loss? 

The problem of production of unobserved per- 
sistent low-mass black-hole X-ray binaries under 
"standard" assuni pt ions was noted before e.g. b y 



Table 1 

Observed short-period low-mass binaries. 



King et all (119961): 



Menou et all (Il999l ): 

Skumanichl(|l972l ) "law" on which 



Ergma and Fedoroval (1998"); 
Ivanova and Kalogera ( 20061 . 



Verbunt and Zwaan 



( 19811 ) model is based, is apparently in conflict 
with observational data on rotation velocities 



in young open cluste rs (|Collier Cameronl 12002 : 
lAndronov et al. . 2003[ ). According to the latter 
study, the time-scale of rotation braking is two or- 
ders of magnitude longer than the one based on 
Skumanich law. Necessity for "weaker" than Ver- 
bunt fc Zwaan AML by MSW was suggested e.g. by 



Hameurv et al.l (|l988l ): llvanova and TaamI ( 2003 ) 



1990) 



for cataclysmic variables and by van Paradijsl 
for SXT. 

In Paper I we suggested, in line with observational 
evidence mentioned above, that magnetic braking 
operates on a much reduced scale (as compared with 



Verbunt and ZwaanI (|l98lh prescription) or that it 



does not operate at all in the semidetached systems 
with black-hole accretors. As a test of this hypoth- 
esis, we computed a population of LMBHB under 
assumption that MSW is not operating once the 
RLOF occurs. 

Since the mechanism that brings the systems to 
the RLOF is the same both in "standard" and in 
"no-MSW" cases, the number of semidetached sys- 
tems is the same in both models (~ 14000), but 
in the latter case, because of a weaker AML, there 
are currently in the Galaxy about 5000 LMBHB 
that evolve to shorter periods and have q > 0.02 
(compared to ~ 3200 in the former case). All these 
systems are transient according to the criterion ([3]) . 
This population is plotted in the lower panel of 
Fig. [2 

Our estimate of the number of SXT apparently ex- 
ceeds the estimates based on observations. Assum- 
ing that their outbursts peak on average at 0.1 icdd 
in X-rays, all of these transients should be seen in 
outburst by the RXTE ASM, regardless of their lo- 



Object 



hr 



Sv 



M© 



1 XTE J1118+480 4.10 K5V-M1V 0.044 - 0.035 6.5 - 7.2 

2 GRO J0422+32 5.09 M2 ± 2V 0.313 - 0.008 3.7 - 5.0 

3 GRS 1009-45 6.84 G5V-M0V 0.159 - 0.125 3.6 - 4.7 

4 XTE J1650-500 7.68 K4 V 0.1 3.4-7.9 

5 A0620-00 7.75 K3V-K7V 0.075 - 0.055 8.7 - 12.9 

6 GS 2000+25 8.28 K3V-K6V 0.053 - 0.035 7.2 - 7.8 

7 XTE J1859-I-226 9.12 G5V-K0V 

8 GRS 1124-68 10.39 K3V-K7V 0.208 - 0.114 6.5 - 8.2 

9 H 1705-25 12.50 K3V-K7V < 0.053 5.6 - 8.3 

cation in the Galaxy. Matching the observed dis- 
covery rate of SXT would require the outbursts to 
be sub-Eddington at maximum and/or recurrence 
times to be ^ 100 years for most sources. 

In the "no-MSW" case, donor-stars in the 
LMBHB have ^ 0.35 at RLOF, while in the 
"standard" case all systems which did not exhaust 
hydrogen in the core to the instant of the RLOF 
evolve to the shorter Porb- This also means that the 
range of initial periods of systems in "no MSW" 
case is smaller: it extends to « 35 hr. 

3.3. Observational parameters of short-period 
LMBHB 

In Table 1 we summarized the main parameters 
of short orbital period SXT that are important from 
evolutionary point of view. The range of the spectral 
types assigned to secondaries is based on published 
data, see for references Paper II; the values of mass 
ratios of components q = M2/Mbh an d the es t imate s 
of black-hole mases are taken from lOrosd ()2003l) : 
(l2004l) . 



Orosz et al. 



3.3.1. Effective temperatures and masses of donors 
The essential information available on the 
companion-stars in SXT is their spectral types. De- 
termination of the latter is a challenging task since 
the emission of the cool star is contaminated by 
the radiation from the accretion disc and the hot 
spot where the accretion stream hits the disc's edge 
( Charles and Coe . 20061 ). Different methods of 5^- 
determination result in a scatter of assigned spectra 
that may amount to several subtypes (Table [T]). In 
the absence of direct determinations of such pa- 
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Fig. 5. Effective temperatures of donors in model population 
for "no-MSW" case (dots). Vertical lines mark the ranges 
of Teff of donor-stars in observed SXT corresponding to the 
ranges of the estimates of their spectral types (Table [!}• 
Systems are annotated according to their numbers in Table[T] 
Large filled circles give T^f{ of donors derived from the fits to 
the synthetic spectra. Dashed lines are evolutionary tracks 
for Mo=1.0 and 1.1 Mq stars with 4 Mq companions . Heavy 
solid lines to the left and right, respectively, show "limiting" 
tracks for (1+12) Mq, Pq = 9.6 hr system in which MSW 
does not operate after RLOF and (1+4) Mq, Pq = 45.6 hr 
system with MSW operating after RLOF (see § 13.3.11 for 
data on tracks and discussion) . Sp — T^g relation used in the 
paper is shown at the right border of the coordinate box. 

rameters of the stars as their effective temperatures 
and masses (for a few exceptions see below), one is 
forced to use Sp — Tcs and Sp — M relations. We 
adopted for the present study Sp — Ten relation by 
Tokuna ga an d Sp — M relation by Schmidt-Kaler 



given m [Cox 



(|2000l ). In Figs. El and m we plot, re- 
spectively, effective temperatures of the donors and 
their masses in model population vs. orbital periods 
of the systems. In the absence of objective criteria 
for discrimination of reliable vs. non-reliable Sp- 
determinations, we compare model data with the 
ranges of Tcff and M2 inferred from the spectra. 

Within uncertainties of the 5p-determinations 
and conversion Sp — Tcff the model satisfactorily re- 
produces Toff and Mi of the donors in the LMBHB 
with Porb < 9 hr. Several systems - XTE 1650-500, 
A0620-00, and GS 2000+25 - seem to be located 
below the "populated" areas of the plots. However, 
we should note the following. In order to avoid 
interpolation between pre-computed tracks that 
evolve upon RLOF continuously to shorter periods 
and tracks that turn to longer periods immedi- 
ately or change the direction of evolution in Porb, 
we restricted ourselves by binaries that evolve to 
shortPorb- For a system with given M\q and M20 
the direction of evolution changes quite abruptly. 



MSW-FGWR pre-contact 
GWR post— contact 
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Fig. 6. Masses of donor-stars in the modeled population 
for "no-MSW" case. Vertical lines show the ranges of M2 
corresponding to the range of spectral types of donors in 
observed SXT (Table[TJ. Dashed and continuous lines are the 
same evolutionary tracks as in Fig. [5] For XTE J1650-500 
(at Porb=7.88 hr) there is only one determination of the 
spectrum - K4V, but the same authors mention that the 
next best fits are G5V and K2III; for this reason we show 
by solid circle the lower limit for M2 in this system. 

over a narrow range of initial APorb ^0.1 day. For 
a given combination of M\q and M20 a "gap" be- 
tween tracks evolving in different directions forms. 
This is seen from Fig. [5] where we plot evolutionary 
tracks for (1+4) M©, Po=1.4 and 1.5 day and for 
(1.1+4) Mq, Po=1.3, 1.4, and 1.45 day systems. 
But since there is a continuity in the initial param- 
eters of the systems, the "gap" has to be filled by 
the systems that start RLOF in the well populated 
area in the Pq — Af2o diagram (Fig. [3]). Thus it is 
at least qualitatively clear that the origin of short- 
period LMBHB may be explained within paradigm 
of strong reduction of magnetic braking in systems 
with donors overflowing Roche lobes. 

In our model we reduced AML by MSW to zero. 
In reality some AML by MSW can be still operat- 
ing. We plot in Figs. [5] and [5] two "limiting" tracks: 
for (1+12) M0, Porb,o = 9.6 hr in which the donor 
overflows Roche lobe almost unevolved and MSW is 
absent after RLOF and for (1+4) Mq, Poib.o = 45.6 
hr in which donor has ~ 10^^ at RLOF and 
MSW continues to operate. Crudely, model popu- 
lations with MSW and without MSW that evolve 
to short Poib have to be located between these two 
limiting curves. There is a contribution of lower and 
higher mass systems of course as we plotted tracks 
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for IM0 donors for simplicity. Adding some AML to 
our model will shift the population to the right, giv- 
ing a better agreement with observations while still 
not producing stable luminous X-ray sources. From 
Fig. 6 in Paper I it is seen that such an adding of 
AML will influence mainly the long-period systems. 

Our computations for the "no-MSW" model show 
that LMBHB evolve to longer periods if Xc ^ 0.35 
at RLOF. For instance a I.IM0 donor with Porbo = 
1.4 (middle curve) spends in the RLOF-state al- 
most 10 Gyr, out of which about 5 Gyr it evolves 
to longer periods. This provides a possibility of ex- 
plaining SXTs with periods > 9-10 hr. 



3.3.2. Mass-transfer rates 

Our model predicts the distribution of SXT over 
mass-transfer rates M that can also be compared to 
M deduced from observations of SXT sources. 

Usually, M for SXT are evaluated by dividing 
the mass accreted during outburst by the recur- 
rence time (see e.g. White and van Paradijsl 1996h . 
Among short-period SXT the recurrence time is 
known only for A0620-00 (about 60 years) and 
4U 1543-47 (about 10 years). For the other systems 
one can only obtain upper limits of AI since only 
one outburst has been seen in the X-ray observa- 
tions epoch (30 years are usually assumed). These 
estimates of M, d erived from the data on outburst 
parameters from IChen et"aL 1 |l997.) and for KV 
U Ma and V406 Vul ba sed o f observation s prese nted 
in IChatv et all (|2003l ) and iHvnes et all (|2002f ) re- 
spectively, are shown as plusses in Figs. H] and [71 

For the "standard" model all above-described M- 
estimates are inconsistent with the model (two up- 
per panels of Fig. |4|). Even if only GWR acts, M ^ 
10~^^ M0 yr~^ cannot be explained if one assumes 
that they represent secular values (lower panel of 
Fig. HI . It would take post-minimum period systems 
several Hubble times to reach Porb ^ 4 hr and have 
a secular M ~ 10~^^ Mp, yr~^. If the extremely low 



M deduced for short period SXTs were due to down- 
ward fluctuations from an intrinsically high secular 
M, one would still have to explain why the bright 
counterparts are not observed fii l3.1|) . Furthermore, 
some systems would then necessarily have higher 
rates than their secular value, hence there should be 
even more persistent systems. 

In Paper I we have shown that in the "standard" 
model systems with Porb ~ 8 — 10 hr may be con- 
sistent with systems where donors overflow Roche 
lobes extremely close to TAMS, but then one needs 



to assume some very special initial distribution of 
binaries over orbital separations which will after sev- 
eral evolutionary stages with mass transfer and mass 
and momentum loss from the system and supernova 
explosion, result in concentration of zero-age black- 
hole and main-sequence-star systems just at the de- 
sired very narrow range of separations. On the other 
hand, if a small amount of AML is added to the al- 
ternative "no-MSW" model, the latter will become 
consistent with recurrence time based estimates of 
M for observed Porb ~ 8 — 10 hr systems, as shows 
comparison of two lower panels of Fig. 3) 

The longest period SXT of the "short-period" 
group, 4U 1543-47 (Porb = 27 hr; it is not shown 
in the Figures) is, in principle, consistent with 
the evolution of "standard" model systems with 
donors overflowing Roche-lobes after formation of 
He-cores. But again, these binaries evolve through 
the Hertzsprung gap extremely fast, and the prob- 
ability ofobserving such systems is very low (see 
iKolbl . Il998l ). In the "no-MSW" model, formation 
of the systems similar to 4U1543-47 is much eas- 
ier explained, since to the range of periods about 
13 hr slowly enough evolve binaries with M20 — 
1 Mq donors which overfilled Roche lobes having 
< 0.35. 



3.3.3. Truncated discs 

It is, however, plausible that the M in short- 
period SXTs are much higher than M derived above. 
These M were obtained treating accretion disc as 
a reservoir which during quiescence is filled up to a 
critical value at which the outburst is triggered and 
the disc emptied. The implicit assumption is that 
the reservoir is not leaky, i.e. that the accretion 
rate at the disc's inner edge is much smaller than 
M. However, the inner truncation of the discs is 
requi red to explain observed quiescent luminosities 
(e.g. Lasota . 19961 Lasota et al. . 1996[ ). The same 
may be true for nq n- quies cer it discs too ( s ee e.g . 
iDone and Gierlihskil . l2006h . lOubus et all (|200ll ) 
showed that the disc instability model can repro- 
duce the observed X-ray light-curves only if discs in 
SXTs are truncated and irradiated. The truncated 
(or "leaky" ) disc paradigm is now commonly used in 
describing the observed timing and spectral proper- 
ties o f X-ray LMBHB (jMcCUntock and Remillard , 
2006L and references therein). 



For systems with non-stationar y quiescent ac cre- 
tion discs M can be estimated as ( LasotaL 2001) 
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M 



D , max 



Mi, 



(5) 



where 



Mo.max - 2.7 X 10'ia-° *^3 



Ml 



-0.38 



-Rd'.iO) (6) 



is the maximum quiescent-disc mass (in g) and e — 
AMD/MD,max " the fraction of the disc's mass lost 
during outburst; a is the kinematic viscosity param- 
eter, Rd,io is the disc's outer radius in 10^° cm and 
Min is the accretion rate at the disc's inner edge. The 
usual, non-leaky disc estimates, neglect Min- How- 
ev er, in SXTs M.n is t he dominant term as argued 
bv lMenou et"all (|l999l ). 

According to the DIM, the disc is in cold thermal 
equilibrium if accretion rate at a ll annuli satisfies 
the inequality I Lasota et al. . 2008[ ) 



M{r) < Ml 



2.64 X 10^\ 



0.01 p2.58 
n 1 -n-in 



-0.85 



g S 



(7) 



where ao.i ~ the viscosity parameter in 0.1. There- 
fore, the mass accretion rate at the truncation ra- 
dius (Min) must be smaller than M ^^■^^{riT^) . Param- 
eterizing the truncation radius as a fraction of the 
circularization radius r\„ = ftVcm (where ft < 0.48, 



see 



Menou et a"l] . ll999[ ) one gets (in Mgyr^^) 



M, 



< 2.5 X 10" 

a/3 



(l + q)'^-'(0.5 - 0.2271og9) 



10.32 



P^'f!-'', (8) 



where period is in days. Inserting in Eq. ([8]) the lower 
limits on q from Table[T]and assuming ft = 0.48, one 
obtains the upper limits to the mass-transfer rates 
shown by arrows in Figs. |4] and [71 One can therefore 
conclude that a substantially reduced MSW AML 
would be consistent with both the stability proper- 
ties and M of SXTs. 

4. Cold, stable systems 

Possible problem of excess of SXT in our model 
mentione d in i) | 3.1 | mav be reso lved in the model sug- 
gested by iMenou et all (|l999l ) (henceforth, MNL). 
MNL pointed out that if the truncation radii were 
slightly larger than estimated to fit observations of 
quiescent SXTs the discs would be globally stable, 
since their accretion rates are lower than the critical 
value given by Eq. ^ corresponding to cold stable 
equilibria. 
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Fig. 7. Mass transfer rates vs. orbital period for the model 
population of semidetached low-mass binaries with assumed 
black hole accretors. Upper panel - model for Verbunt-Zwaan 
assumptions on AML via MSW. Crosses mark systems that 
are stable according to Menou et al. criterion for instability 
of cold discs, while dots mark unstable systems. Middle panel 
shows position of systems with discs stable according to 
Menou et al. criterion in the model in which MSW does not 
operate in semidetached systems. Lower panel — systems with 
unstable discs in model without MSW. Arrows and plusses 
have the same meaning as in Fig. |4] 

Figure [7| shows the distribution of stars in two 
population models with respect to stability crite- 
rion suggested by MNL under assumption the trun- 



cation occurs at r^. 



0.48rc 



In the "standard" 



case with MSW AML aU (~ 2900), but very few 
(~ 60), systems are unstable (Fig. [7] upper panel) 
vs. MNL-criterion. On the other hand, in the model 
without MSW after RLOF there are almost equal 
numbers (~ 2500) of stable and unstable vs. MNL- 
criterion systems (Fig. [71 two lower panels). This 
is because M in the GWR-only scenario are lower 
than in the MSW-f GWR scenario, thereby enabling 
the criterion for a cold disc to be met more easily. 
Truncating the disc at Rdrc would make all tran- 
sients "MNL-stable" for both models. Therefore, it 
is possible that most or even all LMBHB are sec- 
ularly cold/stable systems, erupting from time to 
time due to randomly acting factors like variability 
of truncation radii or mass transfer rate. The pre- 
dicted population size of these systems is too small 
to contribute significantly to the large number of 
faint {Lx ~ 10'^^ erg s^^) hard X-ray sources seen 
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in the d eep Chandra expos ures toward the Galactic 



Center (jMuno et alj . 120031 ) . 



models that are not consistent with observed SXT. 
The issue of acE x A is still controversial (s e e, e.g . 



Tauris and Dewil . l200ll : iPodsiadlowski et alJ . l2003f) 



5. Conclusion 

1. We have modeled populations of short-period 
semidetached close binary systems containing black 
holes and low-mass stars under various assumptions 
about the AML mechanism. 

We have found that the " standard" model (based 



Skumanich stellar rotation braking law 



prescription for AML via magnetic braking) pro- 
duces unobserved persistent LMBHB, while a model 
with pure GWR AML produces only transient sys- 
tems. In the latter model it is possible to reproduce, 
within uncertainty of observations, the number of 
LMBHB in the Galaxy, the effective temperatures 
and masses of the donors in these systems (as in- 
ferred from the spectra of the latter). Mass-transfer 
rates in this model are consistent with the upper 
limits on M , if quiescent accretion discs are trun- 
cated and, hence, leaky. All this suggests that the 
strength of MSW AML should be strongly reduced 
compared to the "standard" model. 

In the model without strong MSW it is unnec- 
essary to invoke systems with the donors that are 
almost at TAMS or already left the latter {Xc < 
0.01) in order to explain observed orbital periods 
of SXT and their mass-exchange rates. More, pop- 
ulation of short-period (Porb ^8 — 9 hr) is formed 
by systems that overflow Roche lobes unevolved or 
slightly evolved {Xc > 0.35), while evolution of sys- 
tems with more evolved donors with masses ~ 1 Mq 
opens possibility to explain the origin of LMBHB 
with longer Porb • 

2. An alternative model suggests that if accre- 
tion discs are truncated in quiescence they could be 
cold and stable. Then truncation of discs at radii 
close to Pcirc can make discs in virtually all LMBHB 
cold and stable. In this case SXT outbursts perhaps 
are not related to DIM and are random events that 
result from "external factors", like episodes of en- 
hanced mass transfer or variability in truncation ra- 
dius. 

3. Population synthesis models presented above 
are obtained under assumption of a high value of 
the product of binding energy parameter of stel- 
lar envelopes and expulsion efficiency of common 
envelopes - acE x A=2. Model with acE x A=0.5 
gives similar results (but a reduced total number 
of systems). Further decrease of acE x A results in 



in the sense that there are no strict criteria for 
defining binding energy of stellar envelopes and 
there is no clear understanding whether sources 
other than gravitational energy may contribute to 
unbinding common envelopes. Apart of common 
envelopes, crucial role in determining the possibility 
of formation of LMBHB is played by stellar wind 
mass-loss in all stages of stellar evolution, since it 
defines parameters of the system at the beginning 
of common envelope stage and immediately before 
supernova explosion that gives birth to black hole. 
If general notions on mass-loss will not be revised, 
our results suggest that existence of short-period 
low-mass X-ray binaries demands high values of 
acE X A parameter. 
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